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Abstract. Previously, several natural integral transforms of the Minkowski question mark 
function F{x) were introduced by the author. Each of them is uniquely characterized by 
certain regularity conditions and the functional equation, thus encoding intrinsic information 
about F{x). One of them - the dyadic period function G{z) - was defined as a Stieltjes transform. 
In this paper we introduce a family of "distributions" Fp{x) for 3? p > 1, such that -^1(2;) is the 
question mark function and F2{x) is a discrete distribution with support on a; = 1. We prove 
that the generating function of moments of Fp{x) satisfies the three term functional equation. 
This has an independent interest, though our main concern is the information it provides about 
F(x). This approach yields the following main result: we prove that the dyadic period function 
is a sum of infinite series of rational functions with rational coefficients. 
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1. Introduction and main result 

The aim of this paper to continue investigations on the moments of Minkowski function, 
begun in [1], [2] and [3]. The function ?(x) ("the question mark function") was introduced by 
Minkowski as an example of a continuous function F : [0, oo) — [0, 1), which maps rationals to 
dyadic rationals, and quadratic irrationals to non-dyadic rationals. For non-negative real x it 
is defined by the expression 

F([ao, ai, a2, ag, •.•]) = 1 - 2-^^" + 2-(''o+"^) - 2-("°+"i+"2) + (1) 

where x = [oq, cii, 02, cis, ...] stands for the representation of a; by a (regular) continued fraction 
|15j . By tradition, this function is more often investigated in the interval [0, 1], and in this case 
it is normalized in order -F(l) = 1, whereas in our case F{1) = |. Accordingly, we make a con- 
vention that 7{x) = 2F{x) for x G [0, 1]. For rational x, the series terminates at the last nonzero 
partial quotient a„ of the continued fraction. This function is continuous, monotone and sin- 
gular [9]. By far not complete overview of the papers written about the Minkowski question 
mark function or closely related topics (Farey tree, enumeration of rationals. Stern's diatomic 
sequence, various 1-dimensional generalizations and generalizations to higher dimensions, sta- 
tistics of denominators and Farey intervals, Hausdorff dimension and analytic properties) can 
be found in [1]. These works include [5], [6], [8], [9], [10], [12], [13] (this is the only paper where 
the moments of a certain singular distribution - a close relative of F{x) - were considered), 
[II], [H], [16], [H], [20], [24], [25] [26], [27], [28], [29], [30], [3l], [33]. The internet page [36] 
contains up-to-date and exhaustive bibliography hst of papers related to Minkowski question 
mark function. 

Recently, in Calkin and Wilf [8] defined a binary tree which is generated by the iteration 

a a a + b 

b a + 6 ' 6 ' 

starting from the root j. The last two authors have greatly publicized this tree, but it was known 
long ago to physicists and mathematicians (alias, Stern-Brocot or Farey tree). Elementary 
considerations show that this tree contains every positive rational number once and only once. 
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each being represented in lowest terms. The first four iterations lead to 




(2) 



It is of utmost importance to note that the nth generation consists of exactly those 2"~^ 
positive rational numbers, whose elements of the continued fraction sum up to n. This fact can 
be easily inherited directly from the definition. First, if rational number | is represented as a 
continued fraction [ ], then the map | — maps | to [oq + 1, Oi..., a^]. Second, the 

map I — > ^ maps | to [0, ai + l,...,aj] in case f < 1, and to [1, oq, ai, a^.] in case f > 1- 
This is an important fact which makes the investigations of rational numbers according to their 
position in the Calkin- Wilf tree highly motivated from the perspective of metric number theory 
and dynamics of continued fractions. 

It is well known that each generation of ([2]) possesses a distribution function and 
Fn{x) converges uniformly to F{x). The function F{ distribution function (in the sense 

of probability theory, which imposes the condition of monotonicity) is uniquely determined by 
the functional equation pQ 

2F(.) = p(--l) + l « -^1. (3) 
\ F(j^) if < i< 1. ^ ' 

This implies F{x) + F{l/x) = 1. The mean value of F{x) has been investigated by several 
authors, and was proved to be 3/2. 

Lastly, and most importantly, let us point out that, surprisingly, there are striking similarities 
and parallels between the results proved in |lj and [2] with Lewis'-Zagier's ([22], [23]) results 
on period functions for Maass wave forms, (see [2] for the explanation of this phenomena). 

Just before formulating the main Theorem of this paper, we provide a short summary of 
previous results proved by the author about certain natural integral transforms of F{x). Let 

oo oo 1 

Ml = jx^dF{x), ^L = j (^^^)^dF(a;) = 2 jx^dF{x). 



Both sequences are of definite number-theoretical significance because 

Mi=lim2^-'^ V [ao,ai,..,a,]^, = lim 2^-" V [0, ai, .., a,]^, 

ao+ai+...+a3=n ai+...+as=n 
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(the summation takes place over rational numbers represented as continued fractions; thus, 
ttj > 1 and ttg > 2). We define the exponential generating functions 



oo 

oo 



L=0 ■ { 





oo 



x). 



One directly verifies that m(t) is an entire function, and that M{t) is meromorphic function 
with simple poles at z = log 2 + 27rm, G Z. Further, we have 



The second identity represents only the symmetry property, given by F{x) + F{l/x) = 1. The 
main result about m(t) is that it is uniquely determined by the regularity condition nx(— t) ^ 
g-Viog2v^^ as t oo, the boundary condition m(0) = 1, and the integral equation 

oo 

m(-s) = (2e" -1) j m'(-t) Jo(2v^) dt, s E R+. (4) 



(Here Jo(*) stands for the Bessel function Jo{z) = ^ JJ^ cos(2;sina;) da;). 

Our primary object of investigations is the generating function of moments. Let G{z) = 

oo 

^ mLZ^~^. This series converges for \z\ < 1, and the functional equation for G{z) (see below) 

L=l 

implies that there exist all derivatives of G{z) at 2; = 1, if we approach this point while remaining 
in the domain 3^2 < 1. Then the integral 

00 1 



(which is Stieltjes transform of F{x)) extends G{z) to the cut plane C\ (1, 00). The generating 
function of moments does not exist due to the factorial growth of M^, but this generating 
function can still be defined in the cut plane C = C \ (0, 00) by dF{x). In fact, this 

integral just equals to G{z + 1). Thus, there exist all higher derivatives of G{z) at 2; = 1, and 
(L-i)! dz^^-i ^('^)L-i ~ ^ 1- The following result was proved in [1]. 



Theorem 1.1. The function G{z), defined initially as a power series, has an analytic contin- 
uation to the cut plane C \ (1, 00) via It satisfies the functional equation 

- + \g(-] +2G{z+l) = G{z), (6) 
z z^ \z/ 
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and also the symmetry property 



g(. + i) = 4g(1 + i)-1, 



(which is a consequence of the main functional equation). Moreover, G{z) — » 0, if z oo and 
the distance from z to a half line [0, oo) tends to infinity. Conversely, the function having these 
properties is unique. 

Accordingly, this result and the specific appearance of the three term functional equation 
justifies the name for G{z) as the dyadic period function. 

We wish to emphasize that the main motivation for previous research was clarification of 
the nature and structure of the moments m^. It was greatly desirable to give these constants 
(emerging as if from geometric chaos) some other expression than the one obtained directly 
from the Farey (or Calkin- Wilf) tree, which could reveal their structure to greater extent. This 
is accomplished in the current work. Thus, the main result can be formulated as follows. 

Theorem 1.2. There exist canonical and explicit sequence of rational functions H„(2;), such 
that for {\z\ < |} U {|z + || < y } one has an absolutely convergent series 



where ^n{z) is polynomial with rational coefficients of degree n — 1. For n > 1 it has the 
following reciprocity property: 



The rational function H„(z) are defined via implicit and rather complicated recurrence ( 127|) 
(see Section 6). The following table gives initial polynomials ^n{z). 



oo 
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n 



1 

--Z 

6 



1 2 2 



2 Q 53 o 

— 2;^ H 2^ 

27 270 



53 

270^ 

4 4 104 n 112 . 224 

— 2 z H 2 z 

81 675 675 2025 



243 



z^ + 



47029 



425250 



1384 
14175' 



787 . 787 



30375 



16 

729' 



1628392 
22325625 



z" + 



272869 
22325625 



2" + 



5392444 
22325625' 



238901 
637875 



z' + 



60750 



477802 
3189375' 



Remark. The constant | can be replaced by any constant less than 1.29^^ (the latter comes 
exactly from Lemma lA.3p . Unfortunately, our method does not allow to prove an absolute 
convergence in the disk \z\ < 1. In fact, apparently the true region of convergence of the series 
in question is the half plane ^z < 1. Take, for example, Zq = | + 4i. Then by ([6]) and symmetry 
property one has 



1 



:G 



G{zo) 
zo ■ 



1 



-G(zo - 1) - — -^g(^— 
1 \ 1 



1 



2(20-2)2 \z0-2J 2(20-1)2 V^o 



:G 



2{zo - 1) 
1 



2(20-2) 2(^0-1)' 



Both arguments under G on the right belong to the unit circle, and thus we can use Taylor 
series for G{z). Using numerical values of m^, obtained via the method described in Appendix 
A.2., we obtain: G(2o) = 0.078083+ + 0.205424+2, with all digits exact. On the other hand, the 
series in Theorem 11.21 for n = 60 gives 

60 

^(-l)"H„(2o) = 0.078090+ + 0.205427+i 

n=0 

Finally, based on the last integral in ([5]), we can calculate G{z) as a Stieltjes integral. If we 
divide the unit interval into = 3560 equal subintervals, and use Riemann-Stieltjes sum, we 
get an approximate value G{zq) ^ 0.078082+ + 0.205424+z. All evaluations match very well. 



Experimental observation 1.3. We conjecture that the series in Theorem converges 
absolutely for 3?2 < 1 . 



With a slight abuse of notation, we will henceforth write /^^ ^H-^o) instead of f^L-i fi 



z=zo 



Corollary 1.4. The moments can he expressed by the convergent series of rational numbers: 



rriL 



lim 2 

n— >oo 



2-n 



00 

[0, ai, a,, a.]^ = — -— $^(-1)"H(^-^)(0) 

ai+a2+...+as=n ^ ' n=0 
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tt{L-1) 
-tin 



(0) 



< 4r, /or every 



The speed of convergence is given by the following estimate: 
M G N. The implied constant depends only on L and M. 

Thus, 7712 = X]^o(~"-'^)"-^n(0) = 0.2909264764+. Regarding the speed, numerical calcula- 
tions show that in fact the convergence is geometric. Theorem 11.21 in case z = 1 gives 

l/2\" 3 



Mi=G(l) = l + + 5^-(- 



n=0 



which we already know (see Corollary 14.51 the above is a Taylor series for Mi( p) in powers of 
p — 2, specialized at po = 1). Geometric convergence would be the consequence of the fact 
that analytic functions rriL^p) extend beyond p = 1 (see below). This is supported by the 
phenomena represented as Experimental observation II. 5[ Meanwhile, we are able to prove only 
the given rate. If we were allowed to use the point z = 1, Theorem 11.21 would give a convergent 
series for the moments as well. This is exactly the same as the series in the Corollary 11.41 
only one needs to use a point z = 1 instead of z = 0. 



Experimental observation 1.5. For L > 1, the series 

1 



(L-1)! 



5^(p-2rH(^-^)(l), M41) 



Mr 



n=0 



has exactly 2 — ^ as a radius of convergence. 

To this account. Proposition 14.31 endorse this phenomena, which is highly supported by nu- 
merical calculations, and which does hold for L = 1. 

The following two tables give starting values for the sequence H^(0). 



77 


h;(o) 


77 


h;(o) 


77 


h;(o) 





1 


5 


7 


10 


8026531718888633 


4 


2 ■ 34 ■ 52 


212 . 39 . 57 . 74 . 11 . 172 


1 





6 


787 


11 


797209536976557079423 


28 ■ 35 ■ 53 


211 . 310 . 58 . 75 . 112 . 173 . 31 


2 


1 


7 


238901 


12 


4198988799919158293319845971 


48 


27 . 36 ■ 54 ■ 7 


214.311.59.76.113.13.174.312 


3 


1 


8 


181993843 


13 


12702956822417247965298252330349561 


~72 




210 . 37 . 55 . 72 


210 . 312 . 510 . 77 . 114 . 132 . 175 . 313 


4 


53 


9 


12965510861 


14 


7226191636013675292833514548603516395499899 


8640 


26 . 38 . 56 . 73 . 17 


216. 313. 511. 78. 115. 133. 176. 314 
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n 


h;(o) 


15 


129337183009042141853748450730581369733226857443915617 
215 . 314 . 512 . 79 . . 134 . 177 . 3^5 . 43 . 127 


16 


31258186275777197041073243752715109842753785598306812028984213251 
218 . 315 . 513 . 710 . . 135 . 178 . 316 . 432 . 1272 


17 


3282520501229639755997762022707321704397776888948469860959830459774414444483 


212 . 316 . 514 . 711 . 118 . 136 . 179 . 317 . 433 . 1273 . 257 



The float values of the last three rational numbers are -0.000025804822076, 0.000018040274062 
and —0.000010917558446 respectively. The alternating sum of the elements in the table is 
En=o(-l)"H;(0) = 0.2909255862+ (where = 17), whereas = 40 gives 0.2909264880+, 
and = 50 gives 0.2909264784+. Note that the manifestation of Fermat and Mersenne primes 
in the denominators at an early stage is not accidental, minding the exact value of the deter- 
minant in Lemma [6. 11 Chapter 6 (see below). Moreover, the prime powers of every odd prime, 
which divides the denominator, increase each time by 1 while passing from H'^(O) to H^^^(O). 
The pattern for the powers of 2 is more complicated. More thorough research of the linear 
map in Lemma 16.11 can thus clarify prime decomposition of denominators; numerators remains 
much more complicated. 



As will be apparent later, the result in Theorem 11.21 is derived from the knowledge of 
p— derivatives of G{ p, z) at p = 2 (see below). On the other hand, since there are two points 
( p = 2 and p = 0) such that all higher p— derivatives of G{ p, z) are rational functions in z, 
it is not completely surprising that the approach through p = also gives convergent series 
for the moments, though in this case we are forced to use Borel summation. At this point, 
the author does not have a strict mathematical proof of this result (since the function G{ p, z) 
is meanwhile deflned only for 3? p > 1), though numerical calculations provide overwhelming 
evidence for its validity. 

Experimental observation 1.6. Define the rational functions (with rational coefficients) 
Qn{z), n>0, by 

Qo(^) = and recurrently by Q„(z) = — - -^Qn-j-i(-l) ■ (^z^ - ^^)- 

Then 

r+l 

^ 00 00 „ / 1 \ f 

= „'™ 2"" E P. «i. = E ( E ^^n^ ■ / *) ^ (7) 

ai+a2 + ...+as=n ^ r=0 n=0 
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Moreover, 

where !^n{z) are polynomials with rational coefficients (Qp integers forp ^ 2) of degree 2n — 2 
with the reciprocity property 



Note the order of summation in the series for rriL, since the reason for introducing exponential 
function is because we use Borel summation. For example, 



"1-2 + 4-8 + 16-32 + ...' 
The following table gives initial results. 



Borel 



oo oo 
r=0 n=0 



:-2) 



r+l 



ee-' dt 



1 

3' 



n 




n 




1 


1 


4 


\{2z^ - 3^5 + Qz* - 3^3 + 6^2 -3^ + 2) 


2 




5 


\{z^ - 2z^ + Az^ - 7z^ + 4z^ - 7z^ + 4^^ - 2z + 1) 


3 




6 


1(2^10 - 5z^ + 12z^ - 202^ + 37z^- 








-20z^ + 37z^ - 20^3 ^ i2z^ - 5z + 2) 



The next table gives Q'„(— 1) = 2(— 1)"^„(— 1) explicitly: these constants appear in the series 
defining the first non-trivial moment m2. Also, since these numbers are p— adic integers for 
p ^ 2, there is a hope for the successful implementation of the idea from the last chapter in 0; 
that is, possibly one can define moments rriL as p— adic rationals as well. 



n q;(-i) 


n QU-l) 


q;(-i) 


q;(-i) 


1 

1 -i 

2 1 

Q _5 

2 

4 25 

5 -16 

6 43 

7 971 
' 8 


8 '-^ 

q 8431 
^ 8 
1 n 50899 
16 

11 -9751 

12 30365 

1 Q 3069719 
32 

1 A 1227099 
4 

1 p; 31719165 
'■^ 32 


1 206836175 
64 

1 7 339942899 
32 

1125752909 

1 n 15014220659 
-"■^ 128 
on 25188552721 
64 

Oi 170016460947 
128 

oo 1153784184807 
256 

OQ 983668214037 
64 


O/i 1685121707817 
32 

op; 92779913448103 
512 

Og 80142274019997 
128 

07 1111839248032133 
512 

00 7740056893342455 

1024 

on 13515970598654393 
512 

20 47354245650630005 

01 665632101181145115 

2048 
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The final table in this section lists float values of the constants 

r+l 

°° Q' (—1 ] r °° 
^r = Yl I ■ / ^"e"* dt, re No, = ^2, 

appearing in Borel summation. 



r 




r 


A- 





0.2327797875 


6 


0.0004701146 


1 


0.0471561089 


7 


0.0004980015 


2 


0.0085133626 


8 


0.0004005270 


3 


0.0005892453 


9 


0.0002722002 


4 


-0.0001872357 


10 


0.0001607897 


5 


0.0002058729 


11 


0.0000812407 



Thus, T.V=o^r = 0.2909400155+ = ma + 0.000013539+. 

This paper is organized as follows. In Section 2, for each p, 1 < p < oo, we introduce 
a generalization of the Farey (Calkin- Wilf) tree, denoted by Qp. This leads to the notion of 
p— continued fractions and p— Minkowski question mark functions Fp{x). Though p— continued 
fractions are of independent interest (one could deflne a transfer operator, to prove an analogue 
of Gauss-Kuzmin-Levy theorem, various metric results and introduce structural constants), we 
confine to the facts which are necessary for our purposes and leave the deeper research for the 
future. In Section 3 we extend these results to the case of complex p, | p — 2| < 1. The crucial 
consequence of these results is the fact that a function X( p, x) (which gives a bijection between 
trees Qi and Qp) is a continuous function in x and an analytic function in p for | p — 2| < 1. In 
Section 4 we introduce exactly the same integral transforms of Fp{x) as was done in a special 
(though most important) case of F{x) = Fi{x). Also, in this section we prove certain relations 
among the moments. In Section 5 we give the proof of the three term functional equation for 
Gp{z) and the integral equation for mp(t). Finally, Theorem 11.21 is proved in Section 6. The 
hierarchy of sections is linear, and all results from previous ones is used in Section 6. Appendix 
A. contains: derivation for the series ([7]); MAPLE codes to compute rational functions H„(z) 
and Qn(^); description of high-precision method to calculate numerical values for the constants 
rriL] auxiliary lemmas for the Section 3. The paper also contains graphs of some p— Minkowski 
question mark functions Fp{x) for real p, and also pictures of locus points of elements of trees 
Qp for certain characteristic values of p. 

2. p-QUESTION MARK FUNCTIONS AND p-CONTINUED FRACTIONS 

In this section we introduce a family of natural generalizations of the Minkowski question 
mark function F{x). Let 1 < p < 2. Consider the following binary tree, which we denote by 
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Qp. We start from the root a; = 1. Further, each element ("root") x of this tree generates two 
"offsprings" by the following rule: 

px X + 1 



We will use the notation Tp(a;) = Up{x) = When p is fixed, we will sometimes discard 
the subscript. Thus, the first four generations lead to 




We refer the reader to the paper [TT], where authors consider a rather similar construction, 
though having a different purpose in mind (see also [6]). Denote by T„( p) the sequence of 
polynomials, appearing as numerators of fractions of this tree. Thus, Ti( p) = 1, T2( p) = p, 
T^Ip) = 2. Directly from the definition of this tree we inherit that 

T2„(p) = pT„(p)forn>l, 
T2„_i( p) = T„_i( p) + p-^T„( p) for n > 2, 

where e = e{n) = 1 if n is a power of two, and e = otherwise. Thus, the definition of 
these polynomials is almost the same as it appeared in [TT] (these polynomials were named 
Stern polynomials by the authors), with the distinction that in [TT] everywhere one has e = 0. 
Naturally, this difference produces different sequence of polynomials. 

There are 2"~^ positive real numbers in each generation of the tree Qp, say a[."\ 1 < k < 2"~^. 
Moreover, they are all contained in the interval [p — 1, -^]- Indeed, this holds for the initial 
root X = 1, and 

1 px 
P - 1 < X < ^ p - 1 < < 1, 

p — 1 X + 1 

1 x + 1 1 
p - 1 < a; < ^ 1 < < 



p- 1 p p- 1 

This also shows that the left offspring is contained in the interval [ p — 1, 1], while the right one 
- in the interval [1, -^^]. The real numbers appearing in this tree have intrinsic relation with 
p— continued fractions algorithm. The definition of the latter is as follows. Let x G ( p — 1, -^)- 
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Consider the following procedure: 

Kx < 
— p-i ' 

p - 1 < a; < 1, 
X = p — 1. 

Then each such x can be uniquely represented as p— continued fraction 

X = [oq, fli, 0,2, as, ....] p, 

where Oj G N for i > 1, and G N U {0}. This notation means that in the course of iterations 
R^lx) we apply T^^(x) exactly ag times, then once X, then we apply exactly ai times, 
then X, and so on. The procedure terminates exactly for those x G ( p — 1, -^zj), which are the 
members of the tree Qp ( " p-rationals" ) . Also, direct inspection shows that if procedure does 
terminate, the last entry > 2. Thus, we have the same ambiguity for the last entry exactly 
as is the case with ordinary continued fractions. At this point it is straightforward to show 
that the nth generation of Qp consists of x = [ag, ai, a^] p such that ^^=0 '^j ~ exactly as 
in the case p = 1 and tree ([2]). 

Now, consider a function Xp(x) with the following property: Xp(x) = x, where x is a rational 
number in the Calkin- Wilf tree ([2]), and x is a corresponding number in the tree ([8]). In other 
words, Xp(x) is simply a bijection between these two trees. First, if x < y, then x <y. Also, all 
positive rationals appear in the tree ([2]) and they are everywhere dense in M+. Moreover, T and 
U both preserve order, and [ p — 1, -^) is a disjoint union of T[ p — 1, -^) and U[p — 1, -^)- 
Now it is obvious that the function Xp(x) can be extended to a continuous monotone increasing 
function 

XpW : R+ ^ [ p - 1, , Xp(oo) = 

Thus, 

Xp([ao, ai, 0^2, as---]) = [^o, oi, 02, ^3---] p- 

As can be seen from the definitions of both trees ((21) and ([H]), this function satisfies functional 
equations 



{T-\x) = px - 1, if 
X(x) = 1, if 
STOP, if 



Xp(x + 1) = 



Xp(x) + 1 
P 

pXp(x) 
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The last one (symmetry property) is a consequence of the first two. We are not aware whether 
this notion of p— continued fractions is new or not. For example, 



[1, 1,1, 1,1,1,. ..]p = Xp(- 



2p L ,,,,,, -^y 2 

= [4,2,1,10,1,1,2,1,5,1,1,2,1,2,1,1,2,1,3,7,4 



2 = [4,l,l,2XT]v2- 

Now fix p, 1 < p < 2. The following proposition follows immediately from the properties of 
F{x). 

Proposition 2.1. There exists a limit distribution of the nth generation of the tree Qp as 
n — >• oo, defined as 

Fp(x) = lim 2-"+^#{A; : 4"^ < x}. 

n— »oo 

This function is continuous, Fp{x) = for x < p — 1, Fp(x) = 1 for x > and it satisfies 
two functional equations: 

f Fp( px - 1) + 1, if 1< X < 
' \ Fp{^), zf p-l<x<l. ^ ^ 

Additionally, 

Fp{x) + ^p(^) = 1 forx>0. 
The explicit expression for Fp{x) is given by 

Fp([ao, ai, a2, ag, ...] p) = 1 - 2-'^« + 2-('^o+"i) - 2-('^«+"i+'^2) + ____ 



We will refer to the last functional equation as the symmetry property. As was said, it is a 
consequence of the other two, though it is convenient to separate it. 

Proof. Indeed, as it is obvious from the observations above, we simply have 

Fp{Xp{x)) = F{x), xe[0,oo). 

Therefore, two functional equations follow from ([3]) and (Q. All the other statements are im- 
mediate and follow from the properties of F{x). □ 

Equally important, consider the binary tree ([H]) for p > 2. In this case analogous proposition 
holds. 
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Proposition 2.2. Let p > 2. Then there exists a limit distribution of the nth generation as 
n ^ oo. Denote it by fp{x) This function is continuous, fp{x) = for x < -^j, fp{x) = 1 for 
X > p — 1, and it satisfies two functional equations: 



2/p(x) 

and 



/p( px - 1) if 1 < X < p - 1, 
/p(x) + /p(^) =l/orx>0. 



Proof. The proof is analogous to the one of Proposition l2.1l only this time we use equivalences 

p - 1 < X < ^ 1 < < P - 1, 

p — 1 X + 1 

1 1 X + 1 

p-l<x< ^ < <p-l. □ 

p-1 p-1 p 

For the sake of uniformity, we introduce Fp{x) = 1 — /p(x) for p > 2. Then Fp{x) satisfies 
exactly the same functional equations ([3]), with a slight difference that -Fp(x) = 1 for x < 
and Fp{x) = for x > p — 1. Consequently, we will not separate these two cases and all our 
subsequent results hold uniformly. To this account it should be noted that, for example, in 
case p > 2 the integral J^ -^-kd-k should be understood as — J^^ ^^d^. Figure 1 gives graphic 
images of typical cases for Fp{x). 



3. Complex case 

After dealing the case of real p, 1 < p < oo, let us consider a tree ([8]) when p G C. For 
our purpose we will concentrate on the case | p — 2| < 1. It should be noted that the method 
which we use allows to extend these result to the case 3?p > 1. The question in determining 
the set in the complex plain where similar results are valid remains open. More importantly, 
the problem to determine all p G C for which there exists an analytic function Gp{z), which 
satisfied (122]) . seems to be much harder and interesting. Here and below [0, oo] stands for a 
compactification of [0, oo). In the sequel, the notion of a function f{z) to be analytic in the 
closed disc I2; — 2| < 1 means that for zq 7^ 1, l^o ^ 2| < 1, this function is analytic in a certain 
small neighborhood of ^o- If 2:0 = 1, this means that there exist all higher derivatives, if one 
approaches the point zq = 1 while remaining in the disc |-2 — 2| < 1. 
In this section we prove the following result. 

Theorem 3.1. There exists a unique function Xp(x) = X( p,x) : {| p — 2| < 1} x {[0, 00]} 

C U {00}, having these properties: 

{i) dC{ p, x) satisfies functional equations 



EXPLICIT SERIES FOR THE DYADIC PERIOD FUNCTION 



15 




p = 1.2, X e [0.2,3] 



p = 3, x e [0.5,2] 




p = 10, a; e [0.1,9] p = 25, a; G [0, 10] 

Figure 1. Functions Fp{x) 



(ii) For fixed p 7^ 1, X{p,x) : [0, 00] — > C is a continuous function, and the image (denote it 

by J^p) is thus a bounded curve; it is contained in the domain {C \ {\z + 1\ < |}; 

{Hi) For every p, | p — 2| < 1, p 7^ 1, m some neighborhood of p there exists the derivative 

p, x), which is a continuous function for x E [0, 00]; 
(iv) There exist all derivatives Sn{x) — P, x) \ p=i : [0, 00) M (the derivatives are taken 

inside | p — 2| < 1). These functions are uniformly continuous for irrational x in any finite 
interval. Moreover, Sn{x) for x > 1, and Sn{x) <SAr 1 for x e (0, 1). 

The curve J^p has a natural fractal structure: it decomposes into two parts, namely 



and ^ "j^ , with a common point z — 1. Additionally, J^p 



As a consequence, ^ J^p for 
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p 7^ 1. Figures 2-4 show the images of J^p for certain characteristic values of p. 

The investigations of the tree Qp deserve a separate paper. I am very grateful to my col- 
leagues Jeffrey Lagarias and Stefano Isola, who sent me various references, also informing about 
the intrinsic relations of this problem with: Julia sets of rational maps of the Riemann sphere; 
iterated function systems; forward limit sets of semigroups; various topics from complex dynam- 
ics and geometry of discrete groups. Thus, the problem is much more subtle and involved than 
it appears to be. This poses a difficult question on the limit set of the semigroup generated by 
transformations Up and Tp, or any other two "conjugate" analytic maps of the Riemann sphere 
(say, two analytic maps A and B are "conjugate", if A{a) = a, B{(3) = (3, A{P) = B{a) for 
some two points a and j3 on the Riemann sphere). Possibly, certain techniques from complex 
dynamics do apply here. As pointed out by Curtis McMuUen, the property of boundedness 
of J^p can be reformulated in a coordinate-free manner. It appears that this curve consists of 
the closure of the attracting fixed points of the elements of the semigroup (Tp,Wp). Then the 
property for the curve being bounded and being bounded away from z = means that it does 
not contain a repelling fixed point of Tp {z = oo) and a repelling fixed point of Up {z = 0). 
It also does not contain neither of the repelling fixed points of the elements of this semigroup. 
Note that 7^(1) = ^2(1) = 1, = ^2(1) — V^- Thus, there exists a small ball D around 

z = 1, such that '72(D) C D, ^2(0) C D, and the last two maps are contractions in D. This 
strict containment is an open condition on p, and thus there exists a neighborhood of p = 2 
such that Theorem 13.11 does hold. I am grateful to Curtis McMullen for this remark: we get 
the result almost for free. Yet, the full result for | p — 2| < 1 is needed. This is not a new 
kind of problem. Some cases of pairs of Mobius transformations were studied. For example, 
the author in [7j deals with the case of a semigroup generated by two maps z sz ± 1, for 
fixed s, \s\ < 1, and investigates a closure of a set of all attracting fixed points. For example, 
for \s\ > 2~^/^ this set is connected. Further development of this problem can be seen in |32j . 
On the other hand, the case of one rational map is rather well understood, and it is treated in 
[4] . Thus, though the machinery of complex dynamics can greatly clarify our understanding of 
the structure of the curve J^p, we will rather employ the techniques from the analytic theory 
of continued fractions. The main source is the monograph by H.S. Wall ^34j- (Lemmas lA.ll 
IA.2l and lA.SI can be found in the Appendix A. 2.) 

Proof of Theorem 13.11 We need the following two results. 

Theorem 3.2. (^ |34] p. 57.) Let v,^, u eN be positive numbers such that 

vi<l, Vu + Vy+i <l, for V >l. (11) 
Suppose given complex numbers e^, z/ G N, such that 

|e,+i|-3?(e,+i)<t;„ z/ > 1. (12) 
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Define the sequence by the recurrence hi = 1, e^+i = -r-r — , u > 1. Then the continued 
fraction 



T= (13) 

1 + ^ 



1 + ^ 



converges if, and only if, (a) some Cy vanishes, or (6) ^ Q for u > 2 and the series Yl'^=i 
diverges. Moreover, if ey{z) : Ki K2 are analytic functions of a complex variable, Ki and 
K2 are compact sets, 01]) and IH^) are satisfied, and the above series diverges uniformly, then 
the continued fraction converges uniformly for all z G Ki . 

Theorem 3.3. (^[34] . p. 60.) If allvy = \, and the conditions (a) and (6) of Theorem \3.2\ hold. 
then -1\<1, ^0. 

For a, 6gN, pGC, |p — 2|<1, define rational functions 

p'' - 1 

^«(P) = -^1 

p + — p 



Since, for fixed p 7^ 1, Wa{ p) ^ p — 1, as a ^ 00, tlien tliere exist two constants ki = ki{ p) 
and k2 = fc2( p), sucli tfiat 

< ki < \Wa{p)\ < k2 < +00, aeN. (14) 

Let X > 1, X = [01,02,03, ...], be an irrational number, Oj G M. Let us consider the continued 
fraction 

1 

^■(p,^) = J^(p, 01,02,. ..) = — ■ (15) 



^ ^ Tg-^^aA P) 
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If x = [ai, a2, aj > 1 is rational, let us define 

1 

jr( p, x) = T{ p, ai, 02, = 



^ai,a2(p) 



1 



1 + Ta^, 

From the definition, this continued fraction obeys the following rule 

1 



p,ai,a2, ...) 



I + Ta^^aii.?) ■ ^i. P, a2,a3---) 
We will now apply Theorem 13.21 to JF( p, ai, 02, as, ...)• Suppose x is irrational. Thus, let 
Cu = Ta^_-^^^a„{ p); > 2. Let US define constants 

fi{a,b)= sup |T,,6(p)| -5R(T,,fe(p)). 

pec,| p-2|<i 

By Lemma [A. 11 fi{a,b) + ^{h^c) < 0.76, a, 6, c G N. Further, from the definition in Theorem 
it follows that 



h2u = WaA?Wa,S P) p"2-l--+-3-a2+ai^ 

= W-^\ p)W^a..^,( P) p'^---'^3+a.-a,_ ^^g) 

It is obvious that the series 'Y^^=i \bv\ diverges. Hence, Theorem 13.21 tells that the continued 
fraction converges, and that for fixed irrational x = [ai,a2,...] > 1, JF( pg, ai, 02, ...) is an 
analytic function in po in some small neighborhood of p. For rational x this is in fact a 
rational function. 

As it is shown in [31], the i^th convergent of the continued fraction (fT3l) (denote it by ^) is 
equal to the z/th convergent (denote it by ^) of the continued fraction 

1 



1 

61 + 



1 

62 + 

1 

&3 + — 

Moreover, since (fTTj) and (fT2|) are satisfied, we have that, for certain positive constant k = 
KhMM) ([3l],p.55-56), 

IQ2.I > A;(l + I62I + I&4I + ...,+|&2.|), 
IQ2.+1I > /c(l + I63I + I&5I + -,+|&2.+i|), (17) 



I QuQu 



+ 1| 
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Now we have 

Proposition 3.4. Fix p G C, | p — 2| < 1, p ^ 1. Let x = [ai,a2, ■■■] > 1 be a real number. 
The function JF( p, x) : [1, oo] — C zs continuous. 

Proof. Fix irrational x > 1. Let 6 > 0, and y > I he sucli that |x — y| < 6. Then there exists 
N such that the first partial quotients of x and y coincide, N = N{6) — > oo as 5 ^ 0. 
Consequently, let the corresponding convergents to JF( p, x) and JF( p, y) be respectively 

.; and 





A2 


A^ 


An+1 


An+2 


Si' 


52 ' 


' Bn' 


Bn+i ' 


Bn+2 


A, 


A2 


An 


A' 


A' 


51 ' 


52' 


' Bn' 


W ' 


B'N+2 



Now, combining IHM . (|T6l) and (|T71) we see that 

IQ2.Q2.+1I > fc'A;?A;2-' X (| p^^ + | p|«3-a2+ai ^ ^ | p^a,^-^-...+a,^a,+a^^ 

p|a.2-ai _j_ I p|a4-a3+a2-ai _j_ _j_ | p|a2iy-----a3+a2-ai 

Denote Ci = k'^klk^^ . Let | p|«2«-i-.-+a3-a2+ai = A^, 1 < £ < z/. The above inequality and the 
arithmetic-harmonic mean inequality give 

ig2.Q2.+ii > ci(Ai + A2 + ... + A,) ■ (I pr^Ar^ + 1 ^rx^' + ... + 1 ^r-K') 

> |p|ci(Ai + A2 + ... + A^)-(Ar^ + A2^ + ... + A;^) > |p|ciz/^ u > 1. 

Analogously we prove that \Q2u-iQ2u\ > I p\c2^'^, > 2. Thus, \QuQu+i\ > cz/^ for certain real 
0, z/ > 2. We see that yield 



Bn\ ^^\QuQu+i\ ^^v^ AT-I 



A 
B 



N 



N 



C 



-1 



This implies p, x) — J^{p,y)\ < In case x is rational we argue in a similar way. In this 
case note that real numbers close to x = [ai, 02, a^] are of the form or [ai, 02, a^, T, ...], 
either [ai, 02, — 1,1,T, ...] for T sufficiently large. The case x = 00 is analogous. This 
establishes the validity of the Proposition. □ 

Eventually, for real number x > 0, x = [oq, ai, a2, ...], let us define 
X( p, [ao,ai,...]) = Wao(p) + 



p-ao 



p-a2 



Wa,{p) + --^ 
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After an equivalence transformation ([34J, p. 19), this can be given an expression 

X( p, [ao, ai, ...]) = Wa,{ p) + p-'^W-i( p) • n p, ai, as, ag, ...)• (18) 

From the very construction, this function satisfies the functional equations ([9]), is continuous 
at a; = 1 and thus is continuous for x G [0, oo]. Obviously, ([9]) determine the values of X( p,x) 
at rational x uniquely, hence a continuous solution to (IHl) is unique. We are left to show that 
the image of the curve J^p is contained outside the circle |-z + 1| < |- This is equivalent to the 
statement that is contained inside the circle \z — p| < But the points on -yrfi are 

exactly the point on the curve J^p with ao = 0. Thus, we need to show that 

I p-iX( p, [0, ai, a2, ...])- 1| = I p-'Vr,-;'^( P, ai, a^, ...) - l\ (19) 

Unfortunately, we cannot apply Theorem I3.3l directly to all p, | p — 2| < 1, since the table above 
Lemma [A . 1 1 shows that /i(l, &) > ^ for infinitely many h. The maximum values /i(l, h) (see the 
definition of this constant) are produced by points p close to x = 2 + e^'^*/'^, or to x- For this 
reason let us introduce 

^J^*{a,b) = sup \Ta,b{ p)| - 3?(T,,fe( p)). 

P6C,| p-2|<l,| p-xl>0.19,| p-xl>0.19 

Then indeed fi*{a,b) < | for all a, 6 G N. Thus, Theorem 13.31 gives |JF( p, ai, a2, ...) — 1| < 1, 
and the statement f|T9l) follows from Lemma IA.3[ In case | p — 2| < 1, | p — x| < 0.19 (or 
I p — xl < 0.19) we use another theorem by Wall ([M], p. 152), which describes the value 
region of a continued fraction (fT3ll . provided elements Cj/ belong to the compact domain in the 
parabohc region \z\ — ^{ze'^'^) < 2/;,cos^|, for certain fixed — vr < < +n, < /i < |. We 
omit the details. This proves part (ii). In a similar fashion we prove part {in). Finally, a 
direct inspection shows that slightly modified proofs remain valid in case p = 1, if we define a 
function to be analytic at p = 1, if it possesses all higher p— derivatives, while remaining inside 
the disc | p - 2| < 1. □ 

Definition 3.5. We define Minkowski p— question mark function Fp{x) : J^p [0,1], by 

Fp(X(p,x)) = F(x), xG[0,cx)]. 

4. Properties of integral transforms of Fp{x) 

For given p, | p — 2| < 1, we define 
p _)_ p"-i _ 2 

Xn = ——n TV' -^n = [Xn, Xn+i] = ^( P, [n, Ti + 1]) for n G Nq. 

P ( P - 1) 

Complex numbers Xn stand for the analogue of non-negative integers on the curve J^p. In 
other words, Xn = W"'(p — 1). We consider as part of the curve J^p contained between 
the points Xn and Xn+i- Thus, Xo = P ~ 1) Xi = li and the sequence Xn is "increasing", in 
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the sense that Xj as a point on a curve J^p is between Xi Xk i < j < k. Moreover, 

oo 

U AU{^}= ^P- 

Proposition 4.1. Let uj{x) : J^p ^ C be a continuous function. Then 

/OO T /* n. 1 

c^(x)dFp(x) = y ^ uf -+ P ~ )dFJx). 



Proof. Indeed, using ( TTOl) we obtain 

oo ^ 



w(x)dFp(x) = ^ / cu(x)dFp(x) = ^ / cu(x)dFp(x)'^ 

oo ^ „ oo _ „ 



n=0 ^ n=0 a 

and this is exactly the statement of the Proposition. □ 
For L, T G No let us introduce 

p"-l \i 



Bl,t{ P) 2J1+1 pTn ( pn.+l 



For example. 



pn 

n=0 ' ' 



Bo,T = TT^f — -, p) 



B4,t{ P) 



2p^-l' "'''^^ (2p^-l)(2p^+i-l)' 

p^(2p^+i + l) 
(2 p^+2 - 1)(2 p^+i - 1)(2 p^ - 1) ' 

pr(4p2T+3^4pT+2^4pT+l^^) 

(2 p^+3 - 1) (2 p^+2 - 1) (2 p^+i - 1) (2 p^ - 1) ' 

pT(2 pT+2 + i)(4 p2T+4 ^ g pT+3 _^ g pT+2 ^ g pT+1 ^ 



(2 p^+4 - 1)(2 p^+3 - 1)(2 p^+2 - 1)(2 p^+1 - 1)(2 p^ - 1) ■ 
As it is easy to see, -Bl,t( p) are rational functions in p for L,T G Nq. Indeed, 

2(S+T) 



^^-^(P) (p_l)L-E pTn2n+l(l pn ) 2( p - 1)^ ' ( , ) E 2" p"(^ 

.y(L\ P- P^^l.t( p) 



(p-1)^ \sj2p'+T-l (2p^+^-l)(2p^+^-i-l)-...-(2p^+i-l)(2p^-l)' 

where ^l,t{ p) are polynomials. This follows from the observation that p = 1 is a root of 
numerator of multiplicity not less than L. 
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As in case p = 1, our main concern are the moments of distributions Fp[x), which are defined 

by 



^0 

1 



= 2 / X\p,x)dF{x) = hm 2^-" V [0, ai, aa, .., a/^., 

Q ai+a2 + ...+as=n 

Ml{p) = J x^dFp(x). 
Thus, if sup^g.y 1^1 = Pp > 1, which is finite for 3?p > 1, p 7^ 1 (see Section 3), then 



Ml(p)<p^. 



Proposition 4.2. T/ie function p) zs analytic in the disc \ p— 2| < 1, including its boundary. 
In particular, if in this disc 



P) •= r = 2^ VvM P - 2) , 



then for any M G N, one has the estimate rj^^i <^ t>~*^ as f —> 00. 

Proof. The function X( p, x) possesses a derivative in pfor3fJp> 1, |p — 2| < 1, and these are 
bounded and continuous functions for x G M+. Therefore rriL^p) has a derivative. For p = 1, 
there exists -^rX( P, x) <C x^~^^, and it is a continuous function for irrational x. Additionally, 
F'{x) = for X G Q+. This proves the analyticity of miip) in the disc | p — 2| < 1. Then an 
estimate for the Taylor coefficients is the standard fact from Fourier analysis. In fact, 

1 





The function mL{2 + e^'^*'^) G C°°(]R), hence the iteration of integration by parts implies the 
needed estimate. □ 



Proposition 4.3. Functions Ml{p) and mL{p) are related via rational functions Bit{p) 'in- 
the following way: 

L 

=0 



L / 
P) = '^'n^s{p)BL-sAP)[ 



Proof. Indeed, this follows from the definitions and Proposition 14. II in case lj{x) = x^. □ 
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Let US introduce, following [I] in case p = 1, the following generating functions: 
mp(t) = X;m^(p)- = 2y e^MFp(x) = y exp(-P^)dFp(x); 



oo 

^ t ^ / ^ dF,(.) ^ / d^.). (20) 

The situation p = 2 is particularly important, since all these functions can be explicitly 
calculated, and it provides the case where all the subsequent results can be checked directly 
and the starting point in proving Theorem ll.2[ Thus, 

t ^ ^ ^ 1 



m2(t) = e% G2{z 



2- z 



By the definition, expressions m£,( p)/ p are Taylor coefficients of G ^{z) at z = 0. Differen- 
tiation of L — 1 times under the integral defining Gp(z), and substitution z = 1 gives 

r 1 °° 
G(^-i)(l) = (^-l)! -<lF,{x) = {L-l)\M,{^)^G,{z + l) = Y,Mdp)z'-\ (21) 

J T n 



^0 



L=0 



with a radius of convergence equal to Pp^. As was proved in [T] and mentioned before, in case 
p = 1 (pi = oo) this must be interpreted that there exist all derivatives at 2; = 1. The next 
Proposition shows how symmetry property refiects in xnp{t). 

Proposition 4.4. One has 

mp(t) = eP*mp(-t). 



Proof. Indeed, 

mp(t) = j exp (^^r^) di^p(a;) = j exp (^pt- ^^7^) dFp(a;) = 

eP*y"exp dFp(x)^=^eP*mp(-t). □ 

This result allows to obtain linear relations among moments p) and the exact value of the 
first (trivial) moment mi(p). 

Corollary 4.5. One has 
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Proof. Indeed, the last propositions implies 

L 

P) = ' 

s=0 

For L = 1 this gives the first statement of the Corollary. Additionally, Proposition 14. 31 for L = 1 
reads as 

and we are done. □ 

5. Three term functional equation 

Theorem 5.1. The function G p{z) can be extended to analytic function in the domain C \ 
( J^p + 1). It satisfies the functional equation 

l + ^GJ^)+2G,iz + l)= pGp(pz), forz^^^^. (22) 

Its consequence is the symmetry property 

G,(.+ l) = -lG,(i + l)-i. 

Moreover, Gp{z) ^ if dist{z, J^p) oo. 

Conversely - the function satisfying this functional equation and regularity property is unique. 
Proof. Let w{x, z) = l_ . Then it is straightforward to check that 



z + 1) = p ■ w{x, pz), 

p 



^ p , 1^ p ^ p^ ^ 

w{ — — -, z+l) = — -w{x, -) - -. 

X + 1 Z'^ z z 

Thus, for |p-2| < 1, p ^ 2, 



2G'p(2 + l) = 2y + l)dFp(x) + 2 j w(a;, ^ + 1) d Fp(a;) 

= / w;(^-,2 + l)dFp(x)+ / + l)dFp(x) 

J X + 1 J p 

.^P 

(In the first integral we used a substitution a; ^ ^). The functional equation holds in case 
p = 2 as well, which can be checked directly. The holomorphicity of G ^{^z) follows exactly as 
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in case p = 1 [Ij. All we need is the first integral in fl20|) and the fact that J" ^ is a closed set. 

As was mentioned, the uniqueness of a function satisfying fl22l) for p = 1 was proved in [1]. 
Thus, the converse implication follows from analytic continuation principle for the function in 
two complex variables ( p, z) (see Lemma [621 below, where the proof in case p = 2 is presented. 
Similar argument works for general p). □ 

Corollary 5.2. Let p 7^ 1, and^ be any closed smooth contour which rounds the curve J^p + l 
once in the positive direction. Then 



Proof. Indeed, this follows from the functional equation fl22l) . as well as from the symmetry 
property. It is enough to take a sufficiently large circle = {\z\ = R} such that + 1 is 
contained in a small neighborhood of 2 = 1, for which ('^"^ + 1) fl ( J^p + 1) = 0. This is 
possible since ^ J^p (see Theorem 13.11) . □ 

We finish with providing an integral equation for mp(t). We indulge in being concise since 
the argument directly generalizes the one used in [1] to prove the integral functional equation 
for m(t) (in our notation, this is mi(t)). 

Proposition 5.3. Let 1 < p < 00 be real. Then the function tnp(t) satisfies the boundary 
condition mp(0) = 1, regularity property mp(— t) <^ e~^*^°^^, and the integral equation 

00 

mp(-s) = y"m'p(-t)(^2eVo(2v^) - Jo(2v^)) dt, seR+. 




For instance, in the case p = 1 this reduces to (jll), and in the case p = 2 this reads as 

00 00 

2e' j e-'Jo{2V2^) dt = 2e'e-^' = + = e"" + j e"* Jo(2v^) dt, 

which is an identity [35] . 

Proof. Indeed, the functional equation for G p{z) in the region < — 1 in terms of m'p(t) reads 
as 

00 



- = / m'{-t)(^—e^ + -e*" - -e^) dt. 
z J '\z+l z z ) 



Now, multiply this by e and integrate over 3fJz = — cr < —1, where s > is real. All the 
remaining steps are exactly the same as in [IJ. □ 
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Remark. If p 7^ 1, the regularity bound is easier than in case p = 1. Take, for example, 
1 < p < 2. Then 



1 

p-1 



1 

p-1 



mp(t)|< j |exp(-^)|dFp(x) < j e'dF^{x)=e'. 



p-1 



p-1 



Thus, Proposition 14.41 gives |mp(— 1)| < e^^ p)*. The same argument shows that for p > 2 we 
have |mp(— 1)| < e"*. 



6. The proof: approach through p = 2 
Let us rewrite the functional equation for G p{z) = G{ p, z) as 

- + p, + 2G{ p,z + l)= pG{ p, p^). (23) 
z z^ \ z / 

Direct induction shows that the following "chain-rule" holds: 

^(pa,p.p.))^E(:)pw'^(-^)-> 

^ i+j=n ^ 

i+j=n-l \ -J / f 

where in the summation it is assumed that i,j > 0. 

Now we will provide rigorous calculations which yield explicit series for G{ p, z) in terms of 
powers of ( p — 2) and certain rational functions. The function G{ p, z) is analytic in {| p — 2| < 
1} X {\A < |}- This follows from Theorem 13.11 and integral representation fl20|) . Thus, for 
(I p — 2| < 1} X {|2;| < |} it has a Taylor expansion 

00 00 

L=l D=0 

Moreover, the function G{2 + e^'^*'^, '^e^-^^v) e C°°(R x R), and it is double-periodic. Thus, 

1 1 



A standard trick from Fourier analysis (using iteration of integration by parts) shows that 
Vv,L <M (4/3)^ ■ (Lv)-^^ for any M G N. Thus, (ES]) holds for ( p, 2) G {| p - 2| < 1} x {j^l < 
3/4}. 

Our idea is a simple one. Indeed, let us look at ( l20l) . This implies the Taylor series for 
niL^p)/ = Y1'^=qVv,l{p — 2y, convergent in the disc | p — 2| < 1. Due to the absolute 



EXPLICIT SERIES FOR THE DYADIC PERIOD FUNCTION 



27 



convergence, the order of summation in fl25l) is not essential. This yields 

oo oo 
i'=0 L=l 

Therefore, let 



n\ (9p' 



-G(p,2;) = H„(z) = Y]r7„,i 

■ p=2 

L=l 



We already know that Ho(^) = gz^- Though rriL^p) are obviously highly transcendental (and 
mysterious) functions, the series for H„(2;) is in fact a rational function in z, and this is the 
main point of our approach. Moreover, we will show that 



where e^„(z) is a polynomial with rational coefficients of degree n — 1 with the reciprocity 
property + 1) = {—ly^z^-'^^nil + -'-)' ^"(0) = 0- We argue by induction on n. First we 
need an auxiliary lemma. 

Let Q[2;]„-i denote the linear space of dimension n of polynomials of degree < n — 1 with 
rational coefficients. Consider a following linear map Cn-i '■ Q[z]n~i — ^ Q[-2]n-i, defined by 



^n-l{P)iz) = Piz + 1) - ^P(2Z) + ^^P(' 



2n+l ^ ' 2"+i \z^ 
Lemma 6.1. det{Cn-i) 0. Accordingly, Cn-i is the isomorphism. 

Remark. Let m = |_|J . Then it can be proved that indeed det(£n-i) = 2^2+^ • 

Proof. Suppose P E ker{Cn~i)- Then a rational function li{z) = ^^^2)"+^ satisfies the three 
term functional equation 



H(2 + 1)-H(22)+h(-)— = 0, z^l. (26) 
Also, ii{z) = 0(1), a.s z —>■ 00. Now the result follows from the next 



Lemma 6.2. Let T{z) be any analytic function in the domain C \ {1}. Then z/ H(z) is a 
solution of the equation 

H(^+1)-H(2;.) + hQ1 = T(z), 
such that ii{z) ^ as z 00, ii{z) is analytic in C \ {2}, then such ii{z) is unique. 
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Proof. All we need is to show that with the imposed diminishing condition, homogeneous 
equation fl26l) admits only the solution H(z) = 0. Indeed, let H(z) be such a solution. Put 
^ ^ 2"z + 1. Thus, 

H(2'^2 + 2) - H(2"+i2 + 2) + — -h( ) = 0. 

^ ^ ^ ' (2"z + 1)2 V2'^z + 1/ 

This is valid for 2; 7^ (since ii{z) is allowed to have a singularity at 2; = 2). Now sum this over 
n > 0. Due to the diminishing assumption, one gets (after additional substitution z ^ z — 2) 



00 ^ 
ll(z) = - V -H 



(2";z - 2"+i + 1)2 V2"z-2"+i + l 

n=0 ^ ^ 

For clarity, put z —z and consider a function H(2;) = H(— 5;). Thus, 

2 



YLiz) = -Y 7 ^^l— 

^ ^ (2"2 + 2"+i - 1)2 V2" 

n=0 ^ ' 



(2"2 + 2"+i - 1)2 V 2"z + 2"+i - 1 

Consider this for z G [0, 2]. As can be easily seen, then all arguments on the right also belong to 
this interval. We want to prove the needed result simply by applying the maximum argument. 
The last identity is still insufficient. For this reason consider its second iteration. This produces 
a series 

H(^) = ^ /2n+m+l^ _^ 2"+''"+2 - 2"Z - 2"+l + 1)2 "^(^™ ° ^"''^ 
n,m=0 ' 

where u!n{z) = 2"2+2"+i-i • said, o c<j„(2;) g [0,2] for z G [0,2]. Since a function H(z) 
is continuous in the interval [0,2], let zq G [0,2] be such that M = \li{zo)\ = sup^g[o,2i |H(-2)|. 
Consider the above expression for z = zq. Thus, 



00 ^ 
M=\U{zo)\< I (2"+™+izo + 2"+™+2 - 2«;zo - 2"+^ + 1)2 ■^('^" ° 



n,m=0 



< 



M V - — = 0.20453+M. 

n,m=0 ^ ' 

This is contradictory unless M = 0. By the principle of analytic continuation, H(z) = 0, and 
this proves the Lemma. □ 



Remark. Direct inspection of the proof reveals that the statement of Lemma still holds with 
a weaker assumption that H(2) is real-analytic function on (— cx),0]. 
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Now, let us differentiate (|23|) n times with respect to p, use f lM|) and afterwards substitute 
p = 2. This gives 

2 (9^' „ /2\ 1 ^1 (9^' /2\ 1 



2H„(^+1)-2H„(2z) + 2hJ-)^. (27) 



zJ z"^ 



We note that this imphes the reciprocity property 

H„(z + 1) = -^H„(i + l), n>l. 

posteriori, this clarifies how the identity F(x) + F{l/x) = 1 reflects in the series for G{z), 
as stated in Theorem 11.21 reciprocity property (non-homogeneous for n = and homogeneous 
for > 1) is reflected in each of the summands separately, whereas the three term functional 
equation heavily depends on inter-relations among H„(z). 

Now, suppose we know all iij{z) for j < n. 
Lemma 6.3. The left hand side of the equation [2^ is of the form 

(z -!)"+!' 

where JJkA G Q[z]„_i. 

Proof. First, as it is clear from the appearance of l.h.s., we need to verify that z does not 
divide a denominator, if l.h.s. is represented as a quotient of two co-prime polynomials. Indeed, 
using the symmetry property in (123|) for the term G( p, ^), we obtain the three term functional 
equation of the form 

1 F 

p — z {p — zy V' p — 2;/ 

Let us perform the same procedure which we followed to arrive at the equation fl27p . Thus, 
differentiation n times with respect to p and substitution p = 2 gives the expression of the 
form 

2 \ 1 



;G(p, — ^) + 2G{ p,z+l)= pG{ p, pz) 

\ D — Z/ 



l.h.s.2 = 2lln{z + 1) - 2H„(2z) - 2H„(- 



2-zJ (2-z)2' 

where lh.s.2 is expressed in terms of iij{z) for j < n. Nevertheless, this time the common 
denominator of l.h.s.2 is of the form {z — l)"'^^{z — 2)"'+^. As a corollary, z does not divide it. 
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Finally, due to the reciprocity property, for n > 1 one has 

This shows that actually l.h.s. = l.h.s.2, and therefore if this is expressed as a quotient of two 
polynomials in lowest terms, the denominator is a power of {z — 1). Finally, it is obvious that 
this exponent is exactly + 1, and one easily verifies that deg ^n{,z) < n — 1. (Possibly, ^n{,z) 
can be divisible by (2; — 1), but this does not have an impact on the result). □ 

Proof of Theorem 11.21 Now, using Lemma 16.11 we inherit that there exists a unique 
polynomial ^n{.z) of degree <n — l such that ^„(-2) = Summarizing, H„(z) = 

(^-2)"+i solves the equation (l27j) . On the other hand. Lemma W2\ implies that the solution of 
f l27|) we obtained is indeed the unique one. This reasoning proves that for|p — 2|<1, |z|<l 
we have the series 



n=0 

This finally establishes the validity of Theorem II. 2[ Note also that each summand satisfies the 
symmetry property. The series converges absolutely for any z, \z\ < 3/4, and if this holds for 
the same does hold for which gives the circle \z + 9/7| < 12/7. □ 

Curiously, one could formally verify that the function defined by this series does indeed 
satisfy (1221) . Indeed, using (I27|) . we get: 



2G{ p, ^ + 1) = 2Ho(2 + 1) + 2 ^( p - 2)"H„(z + 1) 



n=l 

00 / n „ p.j n-l ^ r^j 

2Ho(. + 1) + E( P - 2)" E 71 sijH"-^^^.).^- + J2 -y^H-,-i(2.).^- 

n=l V j=0 j=0 •'' 

i=o •' j=0 

Denote n — j = s. Then interchanging the order of summation for the first term of the sum in 
the brackets, we obtain: 

2 E( P - 2)^^ E 71 9;7H„_,(2z)z^- = 2J2J:(P- 2)^'^> 9;7H.(2z)z^- - 2Ho(2.) = 

n=l j=0 s=0 j=0 

00 

2 J]( p - 2)"H,(2z + ( p - 2)z) - 2Ho(2z) = 2G{ p, p^) - 2Ho(22). 

s=0 
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The same works for the second sum: 

n=l j=0 •'' 

We perform the same interchange of summation for the second and the third summand respec- 
tively. Thus, this yields 



2G( p, ^ + 1) = pG( p, pz) - p, I) + 2Ho(^ + 1) - 2Ho(2^) + ^Ho( 

pG(p, pz) - -^Cfp, -) - -. 

\ z J z 

On the other hand, it is unclear how one can make this argument to work. This would require 
rather detailed investigation of the linear map £n-i and recurrence (l27|l . and this seems to be 
very technical. 

Appendix A. 

A.l. Approach through p = 0. With a slight abuse of notation, we will use the expression 
-^G(0, z) to denote -§^G{ p, z) \ for s G Nq. Though the function G{ p, z) is defined only 
for 3? p > 1, 2; ^ ( + 1), assume that we are able to prove that it is analytic in p in a certain 
wider domain containing a disc | p| < tu, -07 > 0. These are only formal calculations, but they 
unexpectedly yield series ([7]) (see Section 1), and numerical calculations do strongly confirm 
the validity of it. 

Thus, substitution p = into ( l23ll gives G(0, z) = 2{i-~z) ■ Partial differentiation of (l23ll with 
respect to p, and consequent substitution p = gives 

1g(0, 0) + 2^G(0, z+l) = G{0, 0) ^G(0, .) = 

In the same fashion, differentiating the second time, we obtain -j^G(0,z) = ^2(2-1)^^ • 
general, differentiating (1231) n > 1 times with respect to p, using (IMl) . and substituting p = 0, 
we obtain: 



i+j=n—l 



Let 



Then 



;^-^G(0,.) = QJ.). 



n-l rjj 

2Q„(. + 1) = ^-^Q„_,._,(0)(.^- 
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Consequently, we have a recurrent formula to compute rational functions Q(-2). Let Qniz) 
+ Thus, 

Q„W . (i±il(f_il*(£)_ 



where ^„ are polynomials of degree 2n — 2 with the reciprocity property ^n{z) = 
(this is obvious from the recurrence relation which defines Q,n{z))- Moreover, the coefficients 
of are Qp integers for any prime p ^ 2. These calculations yield a following formal result: 



G(p,^)" = "X^p"-Q„(^-l) = X^p' 



„z(^-2)^„(z-l) 



n=0 n=0 ^ ^ 

This produces the "series" for the second and higher moments of the form 

oo 

m2(p)= P^-5^P"Q;(-1). 

n=0 

In particular, inspection of the table in Section 1 (where the initial values for Q^(— 1) are 
listed) shows that this series for p = 1 does not converge. However, the Borel sum is properly 
defined and it converges exactly to the value m2. This gives empirical evidence for the validity 
of d?]). The principles of Borel summation also suggest the mysterious fact that indeed G{ p, z) 
analytically extends to the interval p G [0, 1]. 

Additionally, numerical calculations reveal the following fact: the sequence \/|Qn(— 1)| is 
monotonically increasing (apparently, tends to oo), while ^log |Q^(— 1)| — logn monotonically 
decreases (apparently, tends to — oo). Thus, 

A-<|QU-l)|<(cn)", 

for c = 0.02372 and A = 3.527, n > 150. We do not have enough evidence to conjecture the 
real growth of this sequence. If c = c{n) — > 0, as oo, then the function 

n=0 

is entire, and in case L = 2, result (JTj) is equivalent to the fact that 

A(t)e"*dt = 7712. 







A. 2. Auxiliary lemmas. These lemmas are needed in Section 3. For a,b E N, p G C, 
I p — 2| < 1, define rational functions 



pa+l_pa' "."V^^ a V t, v r (pa_l)(p6_iy 
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Let US define constants 

/i(a,6)= sup |T„,,(p)|-3fJ(T„,,(p)). 

pec,| p-2|<i 

The following table provides some initial values for constants /i(a, b), computed numerically. 



b\a 


1 


2 


3 


4 


5 


6 


1 


0.25000000 


0.01250000 


0.00780868 


0.03343231 


0.05778002 


0.07712952 


2 


0.29846114 


0.03125000 


0.00159908 


0.01212467 


0.02539758 


0.03645721 


3 


0.35999295 


0.05097235 


0.00647895 


0.00676996 


0.01624300 


0.02437494 


4 


0.41433340 


0.07007201 


0.01316542 


0.00500146 


0.01287728 


0.01963810 


o 


yj.'-tooc/yj 1 o t 


08747694 


090fiQ451 


00437959 


01 1 63446 


01 781467 


6 


0.48390408 


0.10255189 


0.02845424 


0.00812804 


0.01125132 


0.01728395 


7 


0.49985799 


0.11503743 


0.03601828 


0.01200557 


0.01120308 


0.01729854 


8 


0.50642035 


0.12494927 


0.04309384 


0.01611126 


0.01125789 


0.01748823 


9 


0.50681483 


0.13248892 


0.04949922 


0.02025219 


0.01132055 


0.01767914 


10 


0.50452450 


0.13796512 


0.05514483 


0.02427779 


0.01136245 


0.01780892 


11 


0.50218322 


0.14173414 


0.06001269 


0.02807992 


0.01138335 


0.01787452 


12 


0.50070286 


0.14415527 


0.06413550 


0.03158969 


0.01139099 


0.01789618 


13 


0.49999979 


0.14555794 


0.06757752 


0.03477145 


0.01139235 


0.01789583 


14 


0.49977304 


0.14622041 


0.07041891 


0.03761547 


0.01139159 


0.01788837 


15 


0.49977361 


0.14636154 


0.07274403 


0.04013040 


0.01139057 


0.01788111 


oo 


0.50000000 


0.12500000 


0.05479177 


0.03097495 


0.01138938 


0.01787406 



Note that there exists lim^^oo Ai(a, 5), and fi{a,b) uniformly in 6, as a — >■ oo. Thus, the 
table above and some standard evaluations give the following 

Lemma A.l. Let a,b,c E N. Then 

li{a, b) + ^l{b, c) < /x(l, 1) + 9) < 0.76. □ 



Lemma A. 2. There exists an absolute constant c > such that for all p e C, 3fJp > 1, and 

all a e N, on has 1-^—^-1 > c. 
' I p-i I 

Proof. Consider a contour, consisting of the segment [1 — iT, 1 + iT], and a semicircle 1 + Te**^, 
— f < < f ■ For sufficiently large T, -^f^ will be large on the semicircle. Moreover, this 
function never vanishes inside or on the contour. Thus, from the maximum-minimum principle, 
its minimal absolute value is obtained on the segment [1 — iT,l + iT]. Thus, let p = ^^iiA^*^, 
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— f < ^/^ < ^. Without loss of generality, let ip >0. Consider the case j- < ip < ^- Then 



p"- 1 



1 2 cos ai/i I 

cos^ ^p 



1 



> 



+ 1 ip'^- 1] 



1 



cos^ V 



- 1 



1 



cosV' 



The function Y{p) is an increasing function in p for p > 1. It is obvious that we may consider 
a case of a sufficiently large. Thus, 



- 1 



> Y 



cos 



iV 



tan 



2 JT_ 

2a 



64a2 



0(1) 
«3 



TT^ I 0(1) 

4a2 a3 



TT^ I 0(1) 

4a2 a3 



r2 
16 



Let now <■?/'< ^. First, consider a function Mog cos(?/'?/') := It is a decreasing 

function for < y < Indeed, this is equivalent to the inequahty 

— tanx ■ X — log cosx < 0, for < x < — . 

The function on the left is itself a decreasing function, with maximum value attained at x = 0. 
Thus, ^(1) > V{a), which means cosatp < cos°"iIj, and this gives 



1 



p-1 



> 



1 



1 



COS 2 -l/l 



- 1 



> 1. □ 



Therefore, Lemma [A. 21 implies that the function p ^W^ ^( p) is uniformly bounded: 

sup I P~^W-\ p)| = Co < +00. 

aeN,| p-2|<l 

This shows the validity of the following Lemma (apart from a numerical bound, which is the 
outcome of computer calculations). 



Lemma A. 3. One has 



sup \p-^W-\p)z-l\ <1.29. □ 

I p-2|<l,aeN,|^-l|<l 



A. 3. Numerical values for the moments. Unfortunately, Corollary IL4I is not very useful 
in finding exact decimal digits of m2. In fact, the vector (mi, m2, ms...) is the solution of 
an (infinite) system of linear equations, which encodes the functional equation ([6]) (see 
Proposition 6). Namely, if we denote cl = Yl'^=i 2^ ~ we have a linear system for nis 

which describes the coefficients rris uniquely: 



TUs = 

L=0 
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Note that this system is not homogeneous {rriQ = 1). We truncate this matrix at sufficiently 
high order to obtain float values. The accuracy of this calculation can be checked on the test 
value mi = 0.5. This approach yields (for the matrix of order 325): 

1712 = 0.2909264764293087363806977627391202900804371021955943665492+, 

1713 = 0.1863897146439631045710466441086804351206556532933915498238+ 
m4 = 0.1269922584074431352028922278802116388411851457617257181016+. 



with all 58 digits exact (note that 3m2 — 2m3 = 0.5). In fact, the truncation of this matrix at 
an order 325 gives rather accurate values for ttil for 1 < L < 125, well in correspondence with 
an asymptotic formula [3j 

rriL = ^47r2 1og2 ■ cq ■ L^/^C^ + 0(C^L-i/^), (28) 

where cq = /g^2'^(l — F{x))dx = 1.030199563382+. So obtained numerical values for higher 
moments tend to deviate from this expression rather quickly. 

Kinney [TB] proved that the Hausdorff dimension of growth points of ?(x) is equal to 

1 



a = i(y"log2(l + x) d?(x))" 





Lagarias [19] gives the following estimates: 0.8746 < a < 0.8749. Tichy and Uitz [33] calculated 
a ^ 0.875. Paradis et al. [26] give the value a ^ 0.874832. We have (note that ?(l-x)+?(a;) = 

1): 



1 1 



A := J log(l + x) d?{x) = j log (l - ^-Y-) ^- (^) + y log 2 d?(x) = 



1 

oo ^ « oo 

- E ^72^ / (1 - + log 2 = - 5^ + log 2. 

L=l Q L=l 

Thus, we are able to present much more precise result: 

^ ^ 0.874716305108211142215152904219159757... 

2A 

with all 36 digits exact. The author of this paper have contacted the authors of [26] inquiring 
about the error bound for the numerical value of a they obtained. It appears that for this 
purpose 10 generations of ([2]) were used. The authors of [26j were very kind in agreeing to 
perform the same calculations with more generations. Thus, if one uses 18 generations, this 
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gives 0.874716 <a< 0.874719. 

Additionally, the constant Cq in fl28l) is given by 

^=0 

This series is fast convergent, and we obtain 

Co = 1.03019956338269462315600411256447867669415885918240... 

A. 4. Rational functions H„ (z). The following is MAPLE code to compute rational functions 
Hn(z)=h[n] and coefficients H'„(0)=alpha[n] for < n < 50. 



> restart; 

> with (Linear Algebra) : 

> U:=50: 

> h[0] :=l/(2-z) : 

> for n from 1 to U do 



> j [n] :=l/2*simplify( 

> add( unapplyCdiff (h[n-j] ,z$j) ,z) (2*z)*2/j !*(z"(j)) ,j=l. .n)+ 

> add( unapplyCdiff (h[n-j-l] ,z$j) ,z) (2*z)*l/j !*(z"(j)) ,j=l. .n-l)+ 

> unapply(h[n-l] ,z) (2*z) ) : 

> k[n] :=simplify((z-l)~(n+l)*(unapply(j [n] ,z) (z)- 

> unapplyCj [n] ,z) (l/z)/z"2)) : 

> M[n, 1] :=Matrix(n,n) :M[n,2] :=Matrix(n,n) : M[n,3] :=Matrix(n,n) : 

> for tx from 1 to n do for ty from tx to n do 

> M[n,l] [ty,tx] :=binomial(n-tx,n-ty) 

> end do : end do : 

> for tx from 1 to n do M[n,2] [tx,tx] :=2"(n-tx) end do: 

> for tx from 1 to n do M[n,3] [tx,n+l-tx] :=2~(tx-l) end do: 

> Y[n] :=M[n,l]-l/2~(n+l)*M[n,2]+(-l)"(n+l)/2"(n+l)*M[n,3] : 

> A[n] :=Matrix(n, 1) : 

> for tt from 1 to n do A[n] [tt, 1] :=coeff (k[n] ,z,n-tt) end do: 

> B [n] : =MatrixMatrixMultiply (Matrixinverse (Y [n] ) , A [n] ) : 

> h[n] :=add(z-(n-s)*B[n] [s,l] (s,l) ,s=l. .n) / (z-2) " (n+1) : 

> end do: 
> 



> for n from to U do alpha[n] :=unapply(dif f (h[n] ,z$l) ,z) (0) end do; 
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It causes no complications to compute h[n] on a standard home computer for < n < 60, 
though the computations heavily increase in difficulty for n > 60. 

A. 5. Rational functions Q„(2;). This program computes Q,n{z) =q[n] and the values 
Q^(-l) =beta[n] for < n < 50. 



> restart; 
>q[0] :=-l/(2*z) ; 
>N:=50: 

>q[l] :=simplify(l/2*unapply(q[0] ,z) (-1) * (l-l/z"2) ) : 

> for n from 1 to N do 

> q[n] :=l/2*simplify( 

> add(unapply(diff (q[n-j-l] ,z$j ) ,z) (-1) /j ! * (z" ( j )-l/z" ( j+2) ) , j=l . . n-l)+ 

> unapply(q[n-l] ,z) (-1) * (l-l/z"2) 

> ): 
end do : 

> for w from to N do beta[w] :=unapply(diff (q[w] ,z$l) ,z) (-1) end do; 



REFERENCES 

[I] G. Alkauskas, The moments of Minkowski question mark function: the dyadic period function (submitted); 
■arXiv: 0801. 0051. 

[2] G. Alkauskas, Generating and zeta functions, structure, spectral and analytic properties of the moments 
of the Minkowski question mark function. Involve (to appear) ; larXiy • 030 1 • 00561 

[3] G. Alkauskas, An asymptotic formula for the moments of Minkowski question mark function in the interval 
[0,1], Lith. Math. J. 48 (4) 2008, 357-367. 

[4] A. F. Beardon, Iteration of rational functions, Springer- Verlag, 1991. 

[5] C. BONANNO, S. Graffi, S. Isola, Spectral analysis of transfer operators associated to Farey fractions, 
Atti Accad. Naz. Lincei CI. Sci. Fis. Mat. Natur. Rend. Lincei (9) Mat. Appl. 19 (1) (2008), 1-23. 

[6] C. BONANNO, S. IsOLA, Orderings of the rationals and dynamical systems, Colloq. Math, (to appear); 
larXiv :0805.217a 

[7] T. BOUSCH, Connexite locale et par chemins holderiens pour les systemes iteres de fonctions; Available at 

[http : / /t opo . math . u-psud . f r/~bousch (1993) (unpublished). 
[8] N. Calkin, H. Wilf, Recounting the rationals, Amer. Math. Monthly 107 (4) (2000), 360-363. 
[9] A. Denjoy, Sur une fonction reelle de Minkowski, J. Math. Pures Appl. 17 (1938), 105-151. 
[10] A. DuSHiSTOVA, N.G. MOSHCHEVITIN, On the derivative of the Minkowski question mark fmiction ?(x); 

arXiv: 0706. 2219 

[II] M.D. ESPOSTI, S. IsOLA, A. Knauf, Generalized Farey trees, transfer operators and phase transitions. 
Comm. Math. Phys. 275 (2) (2007), 297-329. 



38 GIEDRIUS ALKAUSKAS 

fl2 



[13 

[14 

[15 
[16 

[17- 
[18 

[19 

[20" 
[21 
[22 

[23; 

[24 

[25 

[26 

[27" 

[28 
[29 
[30 

[31 

[32 

[33 

[34" 

[35 
[36 



R. GiRGENSOHN, Constructing singular functions via Farey fractions, J. Math. Anal. Appl. 203 (1) (1996), 
127-141. 

P. J. Grabner, p. Kirschenhofer, R. Tichy, Combinatorial and arithmetical properties of linear 
numeration systems, Combinatorica 22 (2) (2002), 245-267. 

M. Kessebohmer, B.O. Stratmann, Fractal analysis for sets of non-differentiability of Minkowski's 
question mark function; J. Number Theory 128 (9) (2008), 2663-2686. 

A. Ya. Khinchin, Continued fractions, The University of Chicago Press, 1964. 

J.R. Kinney, Note on a singular function of Minkowski, Proc. Amer. Math. Soc. 11 (5) (1960), 788-794. 
S. Klavzar, U. Milutinovic, C. Petr, Stern Polynomials. Adv. in Appl. Math. 39 (1) (2007), 86-95. 
J.C. Lagarias, The Farey shift and the Minkowski ?-function, (preprint). 

J. G. Lagarias, Number theory and dynamical systems. The unreasonable effectiveness of number theory 
(Orono, ME, 1991), Proc. Sympos. Appl. Math., 46, Amer. Math. Soc. (1992), 35-72. 

J.C. Lagarias, G.P. Tresser, A walk along the branches of the extended Farey tree, IBM J. Res. 
Develop. 39 (3) May (1995), 788-794. 

M. Lamberger, On a family of singular measures related to Minkowski's ?(a:) function, Indag. Math. 
(N.S.), 17 (1) (2006), 45-63. 

J.B. Lewis, Spaces of holomorphic functions equivalent to the even Maass cusp forms. Invent. Math. 127 
(2) (1997), 271-306. 

J.B. Lewis, D. Zagier, Period functions for Maass wave forms. I, Ann. of Math. (2) 153 (1) (2001), 
191-258. 

H. Okamoto, M. Wunsch, a geometric construction of continuous, strictly increasing singular functions, 
Proc. Japan Acad. Ser. A Math. Set. 83 (7) (2007), 114-118. 

G. Panti, Multidimensional continued fractions and a Minkowski function, Monatsh. Math. 154 (3) (2008), 
247-264. 

J. Paradi's, p. Viader, L. Bibiloni, The derivative of Minkowski's ?(a;) function. J. Math. Anal. Appl. 
253 (1) (2001), 107-125. 

J. Paradi's, P. Viader, L. Bibiloni, A new light on Minkowski's ?(x) function, J. Number Theory 73 
(2) (1998), 212-227. 

G. Ramharter, On Minkowski's singular function, Proc. Amer. Math. Soc, 99 (3) (1987), 596-597. 

S. Reese, Some Fourier-Stieltjes coefficients revisited. Proc. Amer. Math. Soc. 105 (2) (1989), 384-386. 

F. Ryde, On the relation between two Minkowski functions, J. Number Theory, 17 (1) (1983), 47-51. 

R. Salem, On some singular monotonic functions which are strictly increasing. Trans. Amer. Math. Soc. 
53 (3) (1943), 427-439. 

B. SOLOMYAK, On the 'Mandelbrot set' for pairs of linear maps: asymptotic self-similarity, Nonlinearity 
18 (5) (2005), 1927-1943. 

R. F. Tichy, J. Uitz, An extension of Minkowski's singular function, Appl. Math. Lett. 8 (5) (1995), 
39-46. 

H. S. Wall, Analytic theory of continued fractions, D. Van Nostrand Company, Inc., New York, N. Y., 
1948. 

G. N. Watson, A treatise on the theory of Bessel functions, 2nd ed. Cambridge University Press, 1996. 
An exhaustive bibliography on the Minkowski question mark function. Available at 

|http : // www . maths . nott ingheun . ac . uk/ per sonal/pmxga2/ininkowski . htm| 



EXPLICIT SERIES FOR THE DYADIC PERIOD FUNCTION 39 

Vilnius University, The Department of Mathematics and Informatics, Naugarduko 24, Vilnius, 
Lithuania 

Max- Planck- Institut fiir Mathematik, Vivatsgasse 7, 53111 Bonn, Germany 
giedrius . alkauskasQgmail . com 



